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Abstract 

This chapter is devoted to the connection between random matrices and 
maps, i.e graphs drawn on surfaces. We concentrate on the one-matrix model 
and explain how it encodes and allows to solve a map enumeration problem. 

26.1 Introduction 

Maps are fundamental objects in combinatorics and graph theory, originally 
introduced by Tutte in his series of "census" papers |Tut62al ITut62b[ ITut62c|, 
ITut63| . Their connection to random matrix theory was pioneered in the seminal 
paper titled "Planar diagrams" by Brezin, Itzykson, Parisi and Zuber |Bre78j 
building on an original idea of 't Hooft |Hoo74j . In parallel with the idea that 
planar diagrams (i.e maps) form a natural discretization for the random surfaces 
appearing in 2D quantum gravity (see chapter 30), this led to huge develop- 
ments in physics among which some highlights are the solution by Kazakov 
|Kaz86] of the Ising model on dynamical random planar lattices (i.e maps) then 
used as a check for the Knizhnik-Polyakov-Zamolodchikov relations |Kni88j . 

Matrix models, or more precisely matrix integrals, are efficient tools to ad- 
dress map enumeration problems, and complement other techniques which have 
their roots in combinatorics, such as Tutte's original recursive decomposition 
or the bijective approach. Despite the fact that these techniques originate from 
different communities, it is undoubtable that they are intimately connected. 
For instance, it was recognized that the loop equations for matrix models cor- 
respond to Tutte's equations |Eyn06| , while the knowledge of the matrix model 
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result proved instrumental in finding a bijective proof |Bou02bj for a funda- 
mental result of map enumeration |Ben94j . 

In this chapter, we present an introduction to the method of matrix integrals 
for map enumeration. To keep a clear example in mind, we concentrate on the 
problem considered in [Ben94l IBou02b| of counting the number of planar maps 
with a prescribed degree distribution, though we mention a number of possible 
generalizations of the results encountered on the way. Of course many other 
problems can be addressed, see for instance chapter 30 for a review of models of 
statistical physics on discrete random surfaces (i.e maps) which can be solved 
either exactly or asymptotically using matrix integrals. 

This chapter is organized as follows. In section 126.21 we introduce maps 
and related objects. In section 126.31 we discuss the connection itself between 
matrix integrals and maps: we focus on the Hermitian one-matrix model with 
a polynomial potential, and explain how the formal expansion of its free energy 
around a Gaussian point (quadratic potential) can be represented by diagrams 
identifiable with maps. In section 126.41 we show how techniques of random 
matrix theory introduced in previous chapters allow to deduce the solution of 
the map enumeration problem under consideration. We conclude in section [26.5l 
by explaining how to translate the matrix model result into a bijective proof. 

26.2 Maps: definitions 

Colloquially, maps are "graphs drawn on a surface" . The purpose of this section 
is to make this definition more precise. We recall the basic definitions related 
to graphs in subsection 126.2.11 before defining maps as graphs embedded into 
surfaces in subsection 126.2.21 Subsection 126.2.31 is devoted to a coding of maps 
by pairs of permutations, which will be useful for the following section. 

26.2.1 Graphs 

A graph is defined by the data of its vertices, its edges, and their incidence 
relations. Saying that a vertex is incident to an edge simply means that it is 
one of its extremities, hence each edge is incident to two vertices. Actually, 
we allow for the two extremities of an edge to be the same vertex, in which 
case the edge is called a loop. Moreover several edges may have the same 
two extremities, in which case they are called multiple edges. A graph without 
loops and multiple edges is said to be simple. (In other terminologies, graphs are 
always simple and multigraphs refer to the ones containing loops and multiple 
edges.) A graph is connected if it is not possible to split its set of edges into two 
or more non-empty subsets in such a way that no edge is incident to vertices 
from different subsets. The degree of a vertex is the number of edges incident 
to it, counted with multiplicity (i.e a loop is counted twice). 
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(a) 



(b) 



(c) 



(d) 



Figure 26.1: Several cellular embeddings of a same graph into the sphere (drawn 
in the plane by stereographic projection), where half-edges are labelled from 1 
to 10. (a) and (b) differ only by a simple deformation, (c) is obtained by 
choosing another pole for the projection. Thus (a), (b) and (c) represent the 
same planar map encoded by the permutations a = (1 2) (3 4 5) (6 7 8 9) (10) and 
a = (1 3) (2 7) (4 6) (5 9) (8 10). (d) represents a different map, encoded by the 
permutations a' = a o (8 9) and a' = a. Their inequivalence may be checked 
by comparing the face degrees or noting that a o a is not conjugate to a 1 o a'. 

The surfaces we consider here are supposed to be compact, connected, ori- 
entabl^] and without boundary. It is well-known that such surfaces are char- 
acterized, up to homeomorphism, by a unique nonnegative integer called the 
genus, the sphere having genus 0. 

An embedding of a graph into a surface is a function associating each vertex 
with a point of the surface and each edge with a simple open arc, in such a way 
that the images of distinct graph elements (vertices and edges) are disjoint and 
incidence relations are preserved (i.e the extremities of the arc associated with 
an edge are the points associated with its incident vertices). The embedding 
is cellular if the connected components of the complementary of the embedded 
graph are simply connected. In that case, these are called the faces and we 
extend the incidence relations by saying that a face is incident to the vertices 
and edges on its boundary. Each edge is incident to two faces (which might be 
equal) and the degree of a face is the number of edges incident to it (counted 
with multiplicity). It is not difficult to see that the existence of a cellular 
embedding requires the graph to be connected (for otherwise, graph components 
would be separated by non-contractible loops). 

Maps are cellular embeddings of graphs considered up to continuous de- 
formation, i.e up to an homeomorphism of the target surface. Clearly such 
a continuous deformation preserves the incidence relations between vertices, 
edges and faces. It also preserves the genus of the target surface, hence it 
makes sense to speak about the genus of a map. A planar map is a map of 

x It is however possible to extend our definition to non-orientable surfaces. 
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genus 0. In all rigor, this notion should not be confused with that of planar 
graprJl, which is a graph having a least one embedding into the sphere. Figure 
126.11 shows a few cellular embeddings of a same (planar) graph into the sphere 
(drawn in the plane by stereographic projection): the first three are equivalent 
to each other, but not to the fourth (which might immediately be seen by com- 
paring the face degrees). Let us conclude this section by stating the following 
well-known result. 

Theorem 1 (Euler characteristic formula). In a map of genus g, we have 

#{vertices} - #{edges} + #{faces} =2-2g. (26.2.1) 
This quantity is the Euler characteristic \ of the map. 

In the planar case g = 0, the formula is also called the Euler relation. 

26.2.3 Combinatorial maps 

As discussed above, a map contains more data than a graph, and one may 
wonder what is the missing data. It turns out that the embedding into an 
oriented surface amounts to defining a cyclic order on the edges incident to 
a vertexH. More precisely, as we allow for loops, it is better to consider half- 
edges, each of them being incident to exactly one vertex. Let us label the 
half-edges by distinct consecutive integers 1, . . . ,2m (where m is the number 
of edges) in an arbitrary manner. Given a half-edge i, let ot(i) be the other 
half of the same edge, and let a(i) be the half-edge encountered after i when 
turning counterclockwise around its incident vertex (see again Figure 126.11 for 
an example with m = 5). It is easily seen that a and a are permutations of 
{1, . . . ,2m}, which furthermore satisfy the following properties: 

(A) a is an involution without fixed point, 

(B) the subgroup of the permutation group generated by a and a acts tran- 
sitively on {1, . . . , 2m}. 

The latter property simply expresses that the underlying graph is necessarily 
connected. The subgroup generated by a and a is called the cartographic group 
[Zvo95j . 

These permutations fully characterize the map. Actually, a pair of permu- 
tations (a, a) satisfying properties (A) and (B) is called a labelled combinatorial 

2 Unfortunately the literature is often not that careful. To illustrate the distinction let us 
mention that, while the first enumerative results for maps date back to Tutte in the 1960's, 
the mere asymptotic counting of planar graphs (in our present definition) is a fairly recent 
result [Gim05j . 

' ! This observation is generally attributed to Edmonds |Edm60| . For a comprehensive 
graph-theoretical treatment of embeddings into surfaces, we refer to the book of Mohar and 
Thomassen [MohOlj . 
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map and there is a one-to-one correspondence between maps (as denned in the 
previous section) with labelled half-edges and labelled combinatorial maps. In 
this correspondence, the vertices are naturally associated with the cycles of a, 
the edges with the cycles of a and the faces with the cycles of a o a. At this 
stage one may have noticed that vertices and faces play a symmetric role, in- 
deed (cr o a, a) is also a labelled combinatorial map corresponding to the dual 
map. Degrees are given the length of the corresponding cycles, and the Euler 
characteristic is given by 

X(a, a) = c(a) - c(a) + c(a o a) (26.2.2) 

where c(-) denotes the number of cycles. 

Our coding of maps via pairs of permutations depends on an arbitrary 
labelling of the half-edges, which might seem unsatisfactory. We shall identify 
configurations differing by a relabelling, and it easily seen that this amounts to 
identifying (a, a) to all (poaop^ 1 , poaop^ 1 ) where p is an arbitrary permutation 
of {1, ... , 2m}. These equivalence classes are in one-to-one correspondence with 
unlabeled maps. This distinction has some consequences for enumeration: for 
instance the number of labelled maps with m edges is not equal to (2m)! times 
the number of unlabeled maps, because some equivalence classes have fewer 
than (2m)! elements (due to symmetries). By the orbit-stabilizer theorem, the 
number of elements in the class of (cr, a) is (2m)!/r(<7, a) where T(a,a) is the 
number of permutations p such that a = p o a o p~ l and a = p o a o p~ x (such 
p is an automorphism). T(a, a) is often called the "symmetry factor" in the 
literature, though it is seldom properly defined. 

As we shall see in the next section, matrix integrals are naturally related 
to labelled maps. Enumerating unlabeled maps is a harder problem as it re- 
quires classifying their possible symmetries, which is beyond the scope of this 
tex10. The distinction is circumvented when considering rooted maps i.e maps 
with a distinguished half-edge (often represented as a marked oriented edge): 
such maps have no non-trivial automorphism hence the enumeration problem is 
equivalent in the labelled and unlabeled case. Most enumeration results in the 
literature deal with rooted maps. Furthermore maps of large size are "almost 
surely" asymmetric, so the distinction is irrelevant in this context. 

26.3 From matrix integrals to maps 

In this section, we return to random matrices in order to present their con- 
nection with maps. We will concentrate on the so-called one-matrix model 

4 For more on the topic of combinatorial maps and their automorphisms, we refer the 
reader to |Cor75l ICor92] and references therein. [Cor92] also discusses hypermaps, which 
are the natural generalization of combinatorial maps obtained when relaxing constraint (A): 
these are actually bipartite maps in disguise, associated with the two-matrix model of equation 
(|26.3.26|l . 
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(though we will allude to its generalizations): maps appear as diagrams repre- 
senting the expansion of its partition function around a Gaussian point. Our 
goal is to explain this construction in some detail, as this might be also useful 
for the comprehension of other chapters. 

This section is organized as follows. Subsection 126.3.11 provides the defi- 
nitions and the main statement (Theorem [2]) of the topological expansion in 
the one-matrix model. The following subsections are devoted to its derivation 
and generalizations. Subsection 126.3.21 discusses Wick's theorem for Gaussian 
matrix models. Subsection 126.3.31 introduces ah initio the diagrammatic expan- 
sion of the one-matrix model. Subsection 126.3.41 formalizes this computation 
and shows its natural relation with combinatorial maps. Subsection 126.3.51 fi- 
nally presents a few generalizations of the one-matrix model. 

26.3.1 The one-matrix model 

We consider the model of a Hermitian random matrix in a polynomial potential, 
often called simply the one-matrix model, which has been already discussed in 
previous chapters. Different notations and conventions exist, for the purposes 
of this chapter we define its partition function as 

r e NTr(-M 2 /(2t)+V(M)) dM 

where dM is the Lebesgue measure over the space of Hermitian matrices, and 
V stands for a "perturbation" of the form 

oo 

V{x) = ^2—x n . (26.3.2) 

n=l 

Here we shall consider the coefficients v n as formal variables hence, rather than 
a polynomial, V is a formal power series in x and the v n . In this sense, a more 
proper definition of the partition function H^v(i, V) is 

~ N (t,V) = (e NT * v ^) (26.3.3) 

where (•) denotes the expectation value with respect to the Gaussian measure 
proportional to e ~ Nn[tM /( 2t ) dM, acting coefficient-wise on e NTrV ( M ) viewed as 
a formal power series in the v n whose coefficients are polynomials in the matrix 
elements. In other words, the matrix integral in (|26.3.ip must be understood 
in the "formal" sense of chapter 16. We define furthermore the free energy by 

F N (t,V)=logZ N (t,V). (26.3.4) 

The main purpose of this section is to establish the following theorem, which 
is essentially a formalization of ideas present in [Bre781 IHoo74j . 
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Theorem 2 (Topological expansion). The free energy of the one-matrix model 
has the "topological" expansion 

oo 

F N {t,V) = Y,N 2 - 29 F^{t,V) (26.3.5) 

3=0 

where F( 9 \t, V) is equal to the exponential generating function for labelled maps 
of genus g with a weight t per edge and, for all n > 1, a weight v n per vertex 
of degree n. 

Corollary. The quantity 

E<*) (t, V) = 2t dF{9) ^ V) (26.3.6) 

is the generating function for rooted maps (i.e maps with a distinguished half- 
edge) of genus g. Similarly nv n — - corresponds to rooted maps of genus 
g whose root vertex (i.e the vertex incident to the distinguished half-edge) has 
degree n. Maps with several marked edges or vertices are obtained by taking 
multiple derivatives. 

We recall that, from the discussion of section 126.2.31 a labelled map is a 
map whose half-edges are labelled {1, . . . , 2m} where m is the number of edges. 
Hence by exponential generating function for labelled maps we mean 

F®(t,V) = £ j2^y Fi9 ' m) (y) ( 26 - 3 ' 7 ) 



m=0 



where F^< m \V) is the (finite) sum over labelled maps of genus g with m edges 
of the product of vertex weights. If we want to reduce F^ 9 ' m \V) to a sum 
over unlabeled maps instead, then the multiplicity of an individual unlabeled 
map is (2m)!/r, where T is its number of automorphisms. The (2m)! cancels 
the denominator in (|26.3.7p leading to an ordinary generating function, where 
however the weight 1/T has to be kept. It differs from the "true" generating 
function for unlabeled maps where this weight is absent. For rooted maps there 
is no difference since T = 1 for all of them: we simply refer to the generating 
function for rooted maps without specifying between exponential/labelled and 
ordinary/unlabeled. 

F^{t,V) is called the planar free energy. Informally, it is "dominant" in 
the large N limit. Actually, equation (|26.3.5|) makes sense as a sum of formal 
power series (at a given order in t, only a finite number of terms contribute). 

26.3.2 Gaussian model, Wick theorem 

In order to derive theorem we first consider the Gaussian measure 



N 
2wt 



\ N 2 /2 

J e~ NTrM2 /WdM (26.3.8) 
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where dM is the Lebesgue (translation-invariant) measure over the set of Her- 
mitian matrices of size N. It is easily seen that the matrix elements are centered 
jointly Gaussian random variables, with covariance 

{M ij M kl ) = 5 a 5 jk ^. (26.3.9) 

More generally, the expectation value of the product of arbitrarily many matrix 
elements is given via Wick's theorem (generally valid for any Gaussian measure) . 
This classical result can be stated as follows. 

Theorem 3 (Wick's theorem for matrix integrals). The expectation value of 
the product of an arbitrary number of matrix elements is equal to the sum, over 
all possible pairwise matchings of the matrix elements, of the product of pairwise 
covariance. 

For instance, for 4 matrix elements we have 

(M hjl M i2j2 M i3h M i4j4 ) = (M hjl M i2h )(M i3j3 M i4h ) 

+ (M hjl M hh )(M i2]2 M i4ji ) (26.3.10) 
+ {M hjl M i4j4 ){M i232 M i3]3 ). 

Clearly, for an odd number of elements the expectation value vanishes by parity, 
while for an even number 2n of elements the sum involves (2n — 1)!! = (2n — 
l)(2n - 3) • • • 5 • 3 • 1 terms. 

Let us note immediately that these results extend easily to a Gaussian model 
of K random Hermitian matrices . . . , M^ K ^ of same size N, with measure 

{^2^) N ("Y % QaoTrM^M^ dAfW ■ ■ ■ dM™ 

(26.3.11) 

where Q is a K x K real symmetric matrix. The covariance of two matrix 
elements is 

(M^M®) = S a 6 jk ^p^. (26.3.12) 

and, taking the extra index into account, Wick's theorem still applies, for in- 
stance 

(M (ai) M {a2) M {a3) M (ai) ) = (M {ai) M {a2) )(M {a3) M iai) ) 

\ tljl *2J2 1333 »4J4 ' \ lljl 1232 ' V 1333 H34 ' 

+ (mWmWvmWmW) (26.3.13) 

\ lljl 1333 I \ 12J2 «4J4 ' v ' 

+ (M {ai ] M {ai) ) (M {a2 ) M {a3 } ) . 

^ \ Mil U34 I \ 1232 1333 I 
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26.3.3 Diagrammatics of the one-matrix model: a first approach 

We now return to the partition function (|26.3.3p . which is a formal power series 
in the variables v n . If k = (k n ) n >\ denotes a family of nonnegative integers 
with finite support, the coefficient of the monomial v k = Y[^=i 
reads^l 



Inside this expression, each trace Tr M n can be rewritten as a sum of product 
of n elements of M, for instance 

N 

TrM 3 = MijMjkMki. (26.3.15) 

i ,j,k=l 

hence ()26.3.14p may itself be rewritten as the expectation value of a finite linear 
combination of products of elements of M, to be evaluated via Wick's theorem. 



/ x /tt (NTiM n ) kn \ 



3 = k 
- — 



My M i3 M )k M kl (MijMki) 

(a) (b) (c) 

Figure 26.2: Elements constituting the Feynman diagrams for the Hermitian 
one-matrix model: (a) a leg representing a matrix element, (b) a 3-leg vertex 
representing a term in TrM 3 , (c) two paired elements forming an edge. 



We may represent graphically the factors appearing in this decomposition 
as follows (see also Figure l26.2p : 

(a) Each matrix element My is represented as a double line originating from 
a point, forming a leg. The lines are oriented in opposite directions (in- 
coming and outgoing) and "carry" respectively an index i and j. 

(b) Each product of matrix elements appearing in the expansion of Tr M n is 
represented as n legs incident to the same point forming a vertex. The 
legs are cyclically ordered around the vertex so that each incoming line 
is connected to the outgoing line of the consecutive leg, and carries the 
same index. This directly translates the pattern of indices obtained when 



5 A similar expression appears in |Pen88| . albeit with slightly different conventions. 
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writing a trace as product of matrix elements. For instance, for n = 3 the 
decomposition (126.3. 151) yields the vertex shown on Figure [26.2( b). where 
each index i, j, or k takes N possible values. 

(c) Wick's theorem states that the expectation value of a product of matrix 
elements is obtained by matching them pairwise in all possible manners, 
and taking the corresponding product of covariances. A pair of matched 
elements is represented by linking the corresponding legs, forming an edge. 
More precisely, the incoming line of the one leg is connected to the out- 
going line of the other leg, and carries the same index. This translates 
relation (|26.3.9[) : if the connected lines do not carry the same index, then 
the covariance vanishes hence the matching does not contribute to the 
expectation value. 

Globally, the factors appearing in (126.3. 14f) form a collection of k = ^ k n 
vertices, consisting of k n ra-leg vertices for all n > 1. The expectation value is 
obtained by matching the ^ nk n legs, i.e merging them pairwise into edges, in 
all possible manners. Hence, the quantity (126.3. 14h is expressed as a sum over 
all diagrams built out of these vertices and edges, which are sometimes called 
fatgraphs or ribbon graphs in the literature. By the rules discussed above, the 
index lines that are merged together must carry the same index, and they 
form closed oriented cycles. There is clearly a finite number of diagrams, since 
there are (^nA; n )!! ways to merge the legs (in particular, if the number of 
legs is odd, there are no such diagrams, and correspondingly the expectation 
value vanishes). It remains to determine what is the contribution of an indi- 
vidual diagram. Because of relation (|26.3.9|) . each of the m = ^nk n /2 edges 
produces a factor 1/A = t/N. Hence all diagrams will have the same contri- 
bution t m N k ~ m / Yl^Li n kn k n l, taking into account the extra factors present in 
(|26.3.14p . Evaluating this matrix integral amounts to counting the number of 
such diagrams. 




Figure 26.3: The three possible types of diagram appearing in the expansion of 
(TrM 3 • TrM 3 ). 

However, the sum involves many "equivalent" diagrams, i.e diagrams which 
differ only by the choice of line indices or which have the same "shape" . For 
instance, Figure 126.31 displays the three possible types of diagrams obtained 
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when expanding (TrM 3 • TrM 3 ) (in this example, all three diagrams are con- 
nected but this is not true in general). Clearly, we may forget about the line 
indices by counting each index-less diagram with a multiplicity iV , where t is 
the number of cycles of index lines. In our example, we have t = 3 for the first 
two diagrams while I = 1 for the third. Evaluating the shape multiplicity is 
slightly more subtle, and we leave its general discussion to the next subsection. 
It is not difficult to do the computation for the diagrams of Figure 126.31 which 
yields the respective shape multiplicities 9, 3 and 3, and gathering the various 
factors we arrive at 



Generally, as mentioned above, the coefficient (|26.3.14j) will correspond to a 
sum over all (not necessarily connected) diagrams made out of k n vertices with 
n legs for all n. The partition function S^(t,V) is then obtained by summing 
over all k n 's, attaching a weight v n per n-leg vertex, leading to the generating 
function for all diagrams (possibly empty or disconnected). Then F^{t,V) = 
log Hat (t, V) is the generating function for connected diagrams, as it is well- 
known. In Hjv(i, V) as well as in F^it, V), the exponent of N in the contribution 
of a diagram is equal to the number of vertices minus the number of edges plus 
the number of index lines. 



Figure 26.4: The maps corresponding to the diagrams of figure [2531 the first 
two have genus 0, the third genus 1. 



At this stage, it might be rather clear that our (connected) diagrams are 
nothing but maps in disguise (see figure 126.40 . Indeed they are graphs endowed 
with a cyclic order of half-edges (legs) around the vertices, which is a characteri- 
zation of maps as discussed in section 126.21 The cycles of index lines correspond 
to faces, hence the exponent of N in the contribution of a diagram to F^(t, V) 
is equal to the Euler characteristic of the corresponding map. This essentially 
establishes theorem [2j 




(26.3.16) 
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26.3.4 One-matrix model and combinatorial maps 

In this section, we revisit the calculation done above in a more formal manner. 
The purpose is to show that it naturally involves labelled combinatorial maps, 
i.e pairs of permutations, as defined in section [26.2.31 

We start from a "classical" formula for enumerating permutations with pre- 
scribed cycle lengths. Let S p denote the set of permutations of {1, ... ,p} and 
c n (a) denote the number of n-cycles in the permutation a, c(cr) being the total 
number of cycles. Then the numbers of permutations a in S p with prescribed 
values of c n (a) for all n > 1 are encoded into the exponential multivariate 
generating function 



exp 



oo . \ oo 1 oo 



n I <■ — * p. 

\n=l / p=0 a£Sp n=l 



Here (^4 n )n>i is a family of formal variables. Establishing this formula is a 
simple exercise in combinatorics |Fla08j : it simply translates the decomposition 
of a permutation into cycles. We recover Ej^(t,V) of (|26.3.ip on the left-hand 
side by the substitution 

A n = Nv n TrM n (26.3.18) 

and taking the expectation value over the Gaussian random matrix M. By this 
substitution, the a term on the right hand side is, up to a factor independent 
of M, equal to a product of traces which we may rewrite as 

oo p 

jj (TrMT „ M= j2 ri M ^M,)- ( 26 - 3 - 19 ) 

n=l (h,...,i p )e{l,...,N}P q=l 

Wick's theorem and relation (|26.3.9p yield 

\ / t \ p/2 P 

/ V 7 a£l p q=l 

where I p C S p is the set of involutions without fixed point (aka pairwise match- 
ings) of {1, . . . ,p}, which is empty for p odd. We then observe that the product 
on the right hand side of (126.3. 20|) is equal to 1 if the index i q is constant over 
the cycles of the permutation croa, otherwise it is 0. Therefore, when summing 
over all values of . . . , i p ), we find 

/fj(Tr M n ) c "( CT M = (^j^ N c ^ oa \ (26.3.21) 

\n=l / ^ ' a&lp 

Plugging into (|26.3.17p and (|26.3,18p and writing p = 2m = 2c(a) , we arrive at 

oo m oo 
E N (t,V)=^2j— Yl N <-)^(a)+c{aoa) JJ^W. (26.3.22) 

m=0 (ff,a)e5 2m xl 2m n=l 
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We are very close to recognizing a sum over combinatorial maps, with the Euler 
characteristic (|26.2.2p appearing as the exponent of N, but we lack the require- 
ment that a and a generate a transitive subgroup. Again, this is implemented 
by taking the logarithm (in combinatorial terms |Fla08j . 52m x %2m is equal 
to the labelled set construction applied to the class of combinatorial maps - as 
seen by decomposing {1, . . . , 2m} into orbits - and all parameters are inherited, 
hence Hjv(t, V) is the exponential of the generating function for combinatorial 
maps), leading to 

oo m oo 
F N (t, V) = log E N (t, V) = J2 T^yy £ ]] (26.3.23) 

m=0 ^ (o-,o)eA4 m n=l 

where M m is the set of labelled combinatorial maps with m edges. Upon 
regrouping the maps according to their genus, this yields relation (|26.3.7p and 
formally establishes theorem [2J 

26.3.5 Generalization to multi-matrix models 

Let us now briefly discuss multi-matrix models. Informally, we consider a "per- 
turbation" of the multi-matrix Gaussian model (126.3. 11|) by a U (iV)-invariant 
potential N Tr V(M ( - 1 \ , . . .,M^), where V(x 1} x 2 , ... ,x K ) is a "polyno- 
mial" in K non-commutative variables. Actually, V shall be viewed as a formal 
linear combination of monomials in p non-commutative variables x±, x 2 , ■ ■ ■ , xk 
namely 

V(x u x 2 ,...,x K ) = jr / £ V{ai, - n - ,an) x ai ...x an (26.3.24) 

n=0 (ai,...,an)e{l,...,K} n 

where the v(a\, . . . , a n ) are formal variables with the identification v(a 2 , ■ ■ ■ , a n , a{) = 
v(ai, a 2 , . . . , a n ) (since the trace is invariant by cyclic shifts). The partition 
function is defined as the expectation value of e NTrV ( M 1 > M 2 v,A^ (K) ) unc |er 
the Gaussian measure (|26.3.11|) . and the free energy as its logarithm. By a 
simple extension of the arguments of section 126.3.31 (using Wick's theorem for 
multi-matrix integrals), the free energy may be written as a sum over maps, 
whose half-edges now carry one of K "colours" corresponding to the extra in- 
dex a = 1, . . . , K. The formal variable u(oi, . . . , a n ) is the weight for vertices 
around which the half-edge colors are (ax, ■ ■ ■ , a n ) in cyclic order, and is the 
weight per edge whose halves are colored a and b. Furthermore the topological 
expansion (|26.3.5|) still holds. 
In the particular case 

oo K (a) 

V(xx,x 2 ,...,x K ) = Y^T—x n a , (26.3.25) 

n=l a=l 
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all legs incident to a same vertex carry the same colour. The free energy yields 
the generating function for (labelled) maps of a given genus whose vertices are 
colored in K colours, with a weight per vertex with colour a and degree n, 
and a weight Q~\ per edge linking of vertex of colour a to a vertex of colour b. 

Instances of these models appears in several other chapters. Those corre- 
sponding to the form (]26.3.25[) include the chain matrix model in chapter 16, 
the Ising and Potts models in chapter 30. Models of the form (|26.3,24[) but 
not (126.3.25j) include the complex matrix model in chapter 27 (for the diagram- 
matic expansion, M and may be treated as two independent Hermitian 
matrices, which are represented with outgoing and incoming arrows), the 0(n), 
six- vertex and SOS/ADE models in chapter 30. We particularly emphasize the 
two-matrix model with free energy 

J exp f-f Tr M X M 2 + £„ ^ Tr Mf + ^ Tr M») dM 1 dM 2 

log 

/ exp(f TxMiM 2 )dM l dM 2 

(26.3.26) 

which yields generating functions for bipartite maps with a weight per vertex 
depending on degree and colour. It is the most natural generalization of the 
counting problem addressed with the one-matrix model (recovered by setting 

(2) 

Vn = S ni 2/t). 

Let us finally mention that the diagrammatic expansion for real symmetric 
matrices corresponds to maps on unoriented surfaces. 



26.4 The vertex degree distribution of planar maps 

This section is devoted to the enumeration of rooted planar maps with a pre- 
scribed vertex degree distribution, i.e with a given number of vertices of each 
degree. This is equivalent to deriving the generating function for rooted planar 
maps with a weight t per edge and, for all n > 1, a weight v n per vertex of 
degree n. By theorem [2] and its corollary, this in turn amounts to computing the 
derivative with respect to t of the planar free energy of the one-matrix model. 

Let us first state the result, in the form given in [Bou02b] . A different form 
was obtained independently in |Ben94| (without matrices), and a check of their 
agreement can be found in |Bou06] . 

Theorem 4. Let R, S be the formal power series in t and (v n ) n >i satisfying 



R 



S 



(26.4.1) 
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Then the generating function of rooted planar maps with a weight t per edge 
and, for all n, a weight v n per vertex of degree n is given by 

ew = l( R+S 2 -y Vn l y { ^rllt 2){ :r l l ws^^ - 1 

ty f^i P2 j!0"2)!(n-2j + 2)! 

(26.4.2) 

Remark. Clearly, R, S are uniquely determined from the requirement that R = 
S = for t = 0. They have a direct combinatorial interpretation: R is the 
generating function for planar maps with two distinguished vertices of degree 
1 (without weights), S is the generating function for planar maps with one 
distinguished vertex (without weight) of degree 1 and one distinguished face. 

Let us mention that the planar free energy F<°> itself has a more complicated 
expression, involving logarithms. 

This section is organized as follows. In subsection (|26.4.ip we discuss the 
main equation describing the planar limit. In subsection (126. 4. 2p we explain 
how to solve this equation and derive theorem [H Finally in subsection (j26.4.3H 
we present a few instances with explicit counting formulas. 

26.4.1 Saddle-point, loop, Tutte's equations 

Our goal is to "solve" the one-matrix model in the large N limit, i.e extract the 
genus contribution in (|26.3.5 j) . This problem has already been approached 
several times in this book (particularly in chapters 14 and 16), let us recall 
what the "master" equation is. It is an equation is for a quantity called the 
resolvent in the context of matrix integrals. For our purposes, it is nothing but 
a generating function for rooted maps involving an extra variable attached to 
the degree of the root vertex. 

More precisely, we define here the planar resolvent as 

n=0 n=0 n 

where z is a new formal variable and W n is the generating function for rooted 
(i.e with a distinguished half-edge) planar maps whose root vertex (i.e the vertex 
incident to the root) has degree n, with a weight t per edge and, for all n, a 
weight v n per non-root vertex. By convention we set Wo = 1. In comparison 
with the corollary of theorem [2] we do not attach a weight to the root vertex. 
Then, the planar resolvent satisfies the master equation 

W ( z f -{-- V '( z )) W ( z ) + = (26.4.4) 
which is a quadratic equation for I^(z) immediately solved into 

W{z) = \(j- V'{z) ± ^(£-V"(*)) 2 -4P(z)^ . (26.4.5) 
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At this stage P(z) is still an unknown quantity but all derivations of (|26.4.4[) 
show that, unlike W(z), P(z) contains only non-negative powers of z. Actually, 
if V(z) is a polynomial of degree d (i.e we set v n = for n > d), then P{z) is 
a polynomial of degree d — 2. These remarks are instrumental in solving the 
equation, as discussed in the next subsection. Let us first briefly review some 
methods for deriving the master equation. 

Saddle-point approximation. The saddle-point approximation is the orig- 
inal "physical" method used in [Brc78] . It consists in treating the partition 
function (|26.3.1 j) as a "genuine" matrix integral and extracting its analytical 
large N asymptotics. This is done classically by reducing to an integral over 
the eigenvalues, then determining the dominant eigenvalue distribution. See 
chapter 14, sections 1 and 2, for a general discussion of this method. Equation 
(14.2.6) is nothing but equation ()26.4.5p in different notations: the quantities 
denoted by F{z), Q(z) and V'(z) in chapter 14 correspond respectively to W(z), 
P{z) and § - V'(z) here. 

Loop equations. Loop equations correspond to the Schwinger-Dyson equa- 
tions of quantum field theory applied in the context of matrix models |Wad81[ 
Mig83| , see chapter 16 for a general discussion and application in the context 
of the one-matrix model. An interesting feature of loop equations is that they 
provide an easier access to higher genus contributions than the saddle-point 
approximation. However we are here interested in the planar case for which 
they are essentially equivalent. Again, equation (16.4.1) is (j26.4.4|) in different 
notations: the quantities denoted by W^°\x), P^°\z) and V'(z) in chapter 16 
correspond respectively to W(z), P{z) and | — V'{z) here. 
Tutte's recursive decomposition. Tutte's original approach consists in re- 
cursively decomposing rooted maps by "removing" (contracting or deleting) the 
root edge. It translates into an equation determining their generating function, 
upon introducing an extra "catalytic" variable in order to make the decom- 
position bijective. It is now recognized that Tutte's equations are essentially 
equivalent to loop equations [Eyn06| despite their very different origin. 

Let us explain the recursive decomposition in our setting [Tut68, B ou06| . 
We consider a rooted planar map whose root degree (i.e the degree of the root 
vertex) is n, and we decompose it as follows. 

• If the root edge is a loop (i.e connects the root vertex to itself), then it 
naturally "splits" the map into two parts, which may be viewed as two 
rooted planar maps. If there are i half-edges incident to the root vertex 
on one side (excluding those of the loop) , then there are n — 2 — i on the 
other side. These are the respective root degrees of the corresponding 
maps. 

• If the root edge is not a loop, then we contract it (and we may canonically 
pick a new root). If m denotes the degree of the other vertex incident to 
the root edge in the original map, then the root degree of the contracted 
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map is n + m — 2. 



This decomposition is clearly reversible. Taking into account the weights, it 
leads to the equation 

n— 2 oo 

W n = t W l W n _ 2 _ ] + t v m W n+m - 2 (26.4.6) 

i=Q m=l 

valid for all n > 1 with the convention Wo = 1. Equation (|26.4.4p is deduced 
using (|26.4.3|) . in particular P(z) is given by 

oo / oo \ 

P ( z ) = -Y,[ Yl VmW m . 2 . n \z n . (26.4.7) 

n=0 \m=n+2 J 



26.4.2 One-cut solution 

We now turn to the solution of equation (|26.4.4p . In the context of matrix 
models, it gives the "one-cut solution" discussed for instance in chapter 14. Here 
we concentrate on expressing it in combinatorial form. Let us first suppose that 
V(z) (hence P(z)) is a polynomial in z. The one-cut solution is obtained by 
assuming that the polynomial A(z) = (z/t - V'(z)) - AP{z) appearing under 
the square root in (|26.4.5p has exactly two simple zeroes, say in a and b, and 
only double zeroes elsewhere. This leads to 

W(z) = \ (| - V'(z) + G(zW(z-a)(z-b)) (26.4.8) 

where G(z) is a polynomial. This assumption is physically justified in the 
saddle-point picture by saying that the dominant eigenvalue distribution has a 
support made of a single interval [a, b] (corresponding to the cut of W(z)), as 
a perturbation of Wigner's semi-circle distribution. Alternatively, a rigorous 
proof comes via Brown's lemma [Bro65], which translates the fact that W(2) 
hence y A(z) are power series without fractional powers, we refer to |Bou061 
section 10] for details in the current context. y{z — a)(z — b) must be under- 
stood as a Laurent series in 1/z, i.e y/(z — a)(z — b) = z + (lower powers in z). 

Now, it turns out that a, b and G(z) in (|26.4.8|) may be fully determined 
from the mere condition that W(z) = 1/z + (lower powers in z). Indeed, let us 
rewrite (|26.4.8[) as 

'^ V ' [Z) +OW V (26.4.9) 



\/(2-a)(z-6) 2 ^( 2 _ a )( 2 _6) 

Then, we first extract the coefficients of z~ l and z~ 2 on both sides. On the left 
hand side, we obtain respectively and 1 by the above condition. On the right 
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hand side, G(z) does not contribute since it is a polynomial in z. Therefore we 
arrive at 



V'(z) 



y/{z- a)(z - b) 







V{z) 



y/(z- a)(z - b) 



2. 



(26.4.10) 



These equations determine a and b in terms of t and V{z). The coefficients 
may be extracted via a contour integration around z = oo, but a nicer form is 
obtained by performing a change of variable z — > u given by 

R 



u + S + 



u 



(26.4.11) 



also known as Joukowsky's transform. S and R are chosen such that (z—a)(z—b) 
becomes a perfect square namely 

^ 2 

(z-a)(z-b) = ( u 



2 16 
Then, by this change of variable, relations (|26.4,10p yield 



R 

u 



(26.4.12) 



s = t V (u + s + 



R 



R = t + tV'[u + S + 



R 



u 



(26.4.13) 



Upon expanding V'(z) = ^2v n z n , then extracting the respective coefficients 
of vP and u~ l in (u + S + R/u) k via the multinomial formula, we obtain the 
equations (|26.4.ip . 



Extracting further coefficients z 



in (|26.4.9p , we may determine 



step by step the first few W n . In particular, W2 is, up to a factor, the generating 
function E^ for rooted planar maps (without condition on the root vertex), 
since marking a bivalent vertex is tantamount to marking an edge. A slightly 
tedious computation yields 



£(0) 



R + S 2 -V'(u + S+§)\ u _ 3 -2S V ( u + S + §)\ u _ 2 -t 



. (26.4.14) 



which can then be put into the form (|26.4.2p . This establishes theorem [H upon 
noting the restriction that V(z) is a polynomial may now be "lifted": at a 
given order in t, the coefficient of E( ' and the corresponding sum over maps 
both depend on finitely many v n , therefore by suitably truncating V{z) we may 
establish their equality 



26.4.3 Examples 

Tetravalent maps. In the case of a quartic potential V(x) = x 4 /4, the di- 
agrammatic expansion involves only vertices of degree 4, forming 4-regular or 
tetravalent maps. Equations (|26.4.ip and (|26.4.2 j) reduce to 

R-R 3 



E 



(0) 



t 



S = R = t + 3tR\ 



(26.4.15) 
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R is given by a particularly simple quadratic equation solved as 



6t ^ fc + 1 V k 

k=0 v 

where we recognize the celebrated Catalan numbers. Substituting into E^ we 
obtain the series expansion 

^-gaPsji 3 * 4 " <26A17) 

where we identify the number of rooted planar tetravalent maps with 2k edges 
(hence k vertices). This is the same number as that of (general) rooted planar 
maps with k edges [Tut63] . as seen by Tutte's equivalence, and that of rooted 
planar quadrangulations (i.e maps with only faces of degree 4) with k faces, as 
seen by duality. 

Trivalent maps. In the case of a cubic potential V(x) = x 3 /3, the diagram- 
matic expansion involves only vertices of degree 3, forming 3-regular or trivalent 
maps. Equations (|26.4.ip and (126. 4. 2p reduce to 

^(0) = R + S 2 -2R 2 S-t s = t ( S 2 + 2R) R = t + 2tRS (26.4.18) 



and we may eliminate R, yielding a cubic equation for S namely 

(3 _ tsa-tsxi-Bfl 



(hence tS will be a power series in t 3 ). Substituting into £7g , we may use the 
Lagrange inversion formula to compute explicitly its expansion as 



oo 



*f-Er i( ?.>- <*■«<» 

^ (k + 2)\k\\ y ' 

where we recognize the number of rooted planar trivalent maps with 3A: edges 
(hence 2k vertices) |Mul70j . 

Eulerian maps. We now consider the case of a general even potential (y n = 
for n odd). This corresponds to counting maps with vertices of even degree, 
which are called Eulerian (these maps admit a Eulerian path, i.e a path visiting 
each edge exactly once) . A drastic simplification occurs in ()26.4.ip , namely that 
S = and R is given by 

R = t + tJ2^n[ )R n - (26.4.21) 

n=l V n / 

Hence the generating function E^ for rooted planar Eulerian maps depends 
on a single function R satisfying an algebraic equation. This paves the way 
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to an application of the Lagrange inversion formula, allowing to compute the 
general term in the series expansion of E^ namely the number of rooted planar 
Eulerian maps having a prescribed number k n of vertices of degree 2n for all n, 
given by 

2 (En=i"M fr JL ( ™ ~ i) K (26 4 22) 

This formula was first derived combinatorially by Tutte |Tut62cj . 



26.5 From matrix models to bijections 

To conclude this chapter, we move a little away from matrices and explain how 
to rederive the enumeration result of theorem U] through a bijective approach. 
Such an approach consists in counting objects (here maps) by transforming 
them into other objects easier to enumerate. We have already encountered 
several bijections in this chapter: between maps and some pairs of permutations, 
between fatgraphs and maps, in Tutte's recursive decomposition. However they 
do not directly yield an enumeration formula, as a non-bijective step is needed. 

The "easier" objects we shall look for are rooted plane trees (which we 
may view as rooted planar maps with one face). This might not come as a 
surprise ever since the appearance of Catalan numbers at equation (126.4,16j ). 
In general, trees are indeed easy to enumerate by recursive decomposition: 
removing the root cuts the tree into subtrees (forming an ordered sequence due 
to planarity) and this is often immediately translated into an algebraic equation 
for their generating function. Here, we will perform the inverse translation: we 
will construct the trees corresponding to a given equation, namely (126.4.1!) or 
equivalents (126.4.131) . 




Figure 26.5: Left: a 5-tree, together with the canonical matching of black and 
white leaves (dashed lines between arrows). Right: a well-labeled tree. 
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Let us indeed define two classes of rooted trees 1Z and S recursively as 
follows. A 7£-tree (i.e a tree in TV) is either reduced to a "white" leaf, or it 
consists of a root vertex to which are attached a sequence of subtrees that can 
be 7£-trees, 5-trees or single "black" leaves, with the condition that the number 
of such black leaves is equal to the number of 7£-subtrees minus one. A 5-tree 
consists of a root vertex to which are attached a sequence of the same possible 
subtrees, with the condition that the number of black leaves is equal to the 
number of "/^.-subtrees. It is straightforward to check that this is a well-defined 
recursive construction, which translates into the equations (|26.4.ip or ([26.4, 13|) 
for the corresponding generating functions R and S, provided that we attach a 
weight t per vertex or white leaf, and a weight v n per vertex with n — 1 subtrees. 
Figure [26.51 displays a tree in S. 

We may then wonder how such trees are related to maps. A natural idea is 
to "match" the black and white leaves together, creating new edges. Consider 
for instance a 5-tree: it has the same number of black and white leaves, and 
given an orientation there is a "canonical" matching procedure, see again figure 
126.51 This creates a planar map out of a 5-tree, and we further observe that 
the vertex degrees are preserved. It is possible to show that this defines a one- 
to-one correspondence between S and the second class of maps mentioned in 
the remark below theorem HJ see |Bou02bj for details. We proceed similarly 
with 1Z, and then a few more steps allow to establish bijectively theorem HI 
The knowledge of the matrix model solution was instrumental in "guessing" 
the suitable family of trees, which encompasses the one found in [Sch97j based 
on Tutte's formula (126.4.220 . The present construction was further extended to 
bipartite maps corresponding to the two-matrix model (|26.3,26p [Bou02a], and 
maps corresponding to a chain-matrix model [Bou05j . 

The bijective approach has a great virtue. It was indeed realized that it is 
intimately connected with the "geodesic" distance in maps [Bou03] . Actually, 
there is another "dual" family of bijections with so-called well-labeled trees or 
mobiles |Cor8H iMarOl , Bou04[ [Bou07| for which the connection is even more 
apparent. In the simplest instance |Cor81| . a well-labeled tree is a rooted plane 
tree whose vertices carry a positive integer label, in such a way that labels on 
adjacent vertices differ by at most 1 (see again figure [2631 for an example). It 
encodes bijectively a rooted quadrangulation (i.e a map whose faces have degree 
4), a vertex with label I in the tree corresponding to a vertex at distance i from 
the root vertex in the quadrangulation [MarOlj (where the distance is the graph 
distance, i.e the minimal number of consecutive edges connecting two vertices). 
It is easily seen that the generating function for well-labeled trees with root 
label I > 1 satisfies 

R e = t + tRt (Re-i + Re + Ri+i) (26.5.1) 

together with the boundary condition Rq = 0, t being a weight per edge or 
vertex. Through the bijection, Ri yields the generating function for quadran- 
gulations with two marked points at distance at most £, related to the so-called 
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two-point function [Amb95j , Equation (|26.5.ip is nothing but a refinement of 
the third relation of (126.4. 15() (and we have E^ = R\/t). Remarkably, it 
has an explicit solution |Bou031 section 4.1]. It furthermore looks surprisingly 
analogous to (yet different from) the "first string equation" (14.2.14). This 
still mysterious analogy is much more general, as one may refine equations 
(|26.4.ip into discrete recurrence equations involving the distance, similar to the 
string equations for the one-matrix model, and having again explicit solutions 
|Bou03llDiF05] . 

The correspondence between maps and trees has sparked an active field 
of research between physics, combinatorics and probability theory, devoted to 
the study of the geometry of large random maps, see for instance jMie09| and 
references therein. 
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